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Recently, extended phase space thermodynamics of Lovelock AdS black holes has been of great 
interest. To provide insight from a different perspective and gain a unified phase transition picture, 
non-extended phase space thermodynamics of (n+ l)-dimensional charged topological Lovelock AdS 
black holes is investigated detailedly in the grand canonical ensemble. Specifically, the specific heat 
at constant electric potential is calculated and phase transition in the grand canonical ensemble is 
discussed. To probe the impact of the various parameters, we utilize the control variate method 
and solve the phase transition condition equation numerically for the case k = 1,-1. There are two 
critical points for the case n — 6, k = 1 while there is only one for other cases. For k = 0, there 
exists no phase transition point. To figure out the nature of phase transition in the grand canonical 
ensemble, we carry out an analytic check of the analog form of Ehrenfest equations proposed by 
Banerjee et al. It is shown that Lovelock AdS black holes in the grand canonical ensemble undergo a 
second order phase transition. To examine the phase structure in the grand canonical ensemble, we 
utilize the thermodynamic geometry method and calculate both the Weinhold metric and Ruppeiner 
metric. It is shown that for both analytic and graphical results that the divergence structure of the 
Ruppeiner scalar curvature coincides with that of the specific heat. Our research provides one more 
example that Ruppeiner metric serves as a wonderful tool to probe the phase structures of black 
holes. 

PACS numbers: 04.70.Dy, 04.70.-s 


I. INTRODUCTION 


In our recent paper [lj , P — V criticality of topological AdS black holes in Lovelock-Born-Infeld gravity has been 
investigated in the extended phase space and some unique phenomena have been found. It was shown that P — V 
criticality exists not only for the spherical topology but also for k = — 1. This result is really intriguing that it 
has attracted further investigation 2|- Q. On the other hand, it would also be interesting to probe this issue in 
the non-extended phase space to search for some more unique characteristics due to Lovelock gravity. Lovelock 
gravity Q is a particular higher curvature gravity theory which successfully solves the problem of fourth order field 
equations and ghost. In Lovelock gravity, the field equation is only second order and the quantization is free of 
ghosts f?]. Both the black holes and their thermodynamics in Lovelock gravity @-[36} have attracted considerable 
attention. Concerning the thermodynamics of Lovelock black holes in the non-extended phase space, some efforts 
have been made. Topological black hole solutions in Lovelock-Born-Infeld gravity were proposed in Ref. @. Both the 
thermodynamics of asymptotically AdS rotating black branes with flat horizon and asymptotically flat black holes for 
k = 1 were detailedly investigated there. For charged topological AdS black holes, Ref. @ presented the expression of 
the temperature. Ref. [23|,[2J] further studied their entropy and specific heat at constant charge. Ref. [35} studied their 
specific heat and critical exponents in the canonical ensemble. The above research was carried out in the canonical 
ensemble, leaving the grand canonical ensemble unexplored. In this paper, we would like to complete the phase 
transition research of Lovelock charged topological AdS black holes in the grand canonical ensemble. 

In traditional thermodynamics, one can utilize Clausius-Clapeyron-Ehrenfest’s equations to probe the nature of 
phase transitions. The Clausius-Clapeyron equation holds for a first order phase transition while Ehrenfest’s equations 
are satisfied for a second order phase transition. Recently, Banerjee et al. introduced a novel Ehrenfest scheme to 
investigate phase transitions of black holes in the grand canonical ensemble [37j|-[42l|. We utilized this scheme in the 
case of charged topological black hole in Horava-Lifshitz gravity [Hf and also generalized it to the extended phase 
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space 0-0. Ref. [46] further generalized it to the full phase space. The original Ehrenfest equations in traditional 
thermodynamics were utilized in the extended space of black holes in Lovelock-Born-Infeld gravity to study the nature 
of phase transition at the critical point [ 2 }. However, in this paper, we would like to utilize the analog form of Ehrenfest 
scheme proposed by Banerjee et al. to investigate the nature of phase transition points of Lovelock AdS black holes 
in the grand canonical ensemble. 

Different from the traditional thermodynamic method, thermodynamic geometry has served as an alternative way 
to investigate phase transitions of black holes. The well-known examples are Weinhold geometry [47] and Ruppeiner 
geometry (48[. Weinhold defined metric structure in the energy representation as g W = didjM(U, N a ). Here, U is the 
internal energy U while N a represents the extensive thermodynamic variables. Ruppeiner pro pos ed metric structure 
as the Hessian of the entropy. Namely, gF = —didjS(U, N a ). Recently, Quevedo et al. ^ proposed another 
thermodynamic geometry method named as geometrothermodynamics (GTD). For its profound physical meaning, 
Ruppeiner’s metric has been applied to investigate various thermodynamic systems including black holes. For a nice 
review of Ruppeiner geometry, see Ref. [50]. For recent papers, see Ref. & [62| . However, thermodynamic geometry 
of Lovelock AdS black holes in the grand canonical ensemble is still absent in literature. In this paper, we would 
like to explore the Ruppeiner geometry of (n+ l)-dimensional topological AdS black holes in Lovelock gravity in the 
grand canonical ensemble . 

In Sec. m the thermodynamics of charged topological AdS black holes in Lovelock-Born-Infeld gravity will be 
briefly reviewed and the phase transition in the grand canonical ensemble will be investigated in detail. To probe 
the nature of phase transition in the grand canonical ensemble, an analytic check of the analog form of Ehrenfest 
equations will be carried out in Sec. m In Sec. hyi thermodynamic geometry will be studied to examine the phase 
structure of topological AdS black holes. Concluding remarks will be presented in Sec. ]V] 


II. PHASE TRANSITION IN THE GRAND CANONICAL ENSEMBLE 


The action of third order Lovelock-Born-Infeld gravity reads Q 
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/3, a 2 and 03 are Born-Infeld parameter, the second and third order Lovelock coefficients respectively. L(F) denotes 
the Born-Infeld Lagrangian with T4, = dpA v — d v Ap, where Ap is electromagnetic vector. The (n + l)-dimensional 
static solution was derived in Ref. [8] as 
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k and m are parameters related to the curvature of hypersurface and the mass respectively, dfl 2 denotes the line ele¬ 
ment of (n— l)-dimensional hypersurface with constant curvature (n—l)(n — 2)k and F{jg) denotes the hypergeometric 
function as follow 
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where 
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Note that the above solution was derived for the special case that 
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When /? —> 00 , the Born-Infeld Lagrangian reduces to the Maxwell form and the solutions become Lovelock AdS 
black holes. To concentrate on the effects of the third order Lovelock gravity, we will mainly consider Lovelock AdS 
black holes in this paper. 

When /3 —>■ 00 , one can obtain 


, . 3am 6 aA 3 aq 2 

9{r) 1 + + n(n — 1 ) " 


(13) 


The horizon radius r + can be derived from the largest root of the equation f(r) = 0. One can express m in the 
function of r + as 
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Then the mass of (n + l)-dimensional topological AdS black holes can be derived as 
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where E* denotes the volume of the (n — l)-dimensional hypersurface. The Hawking temperature has been derived 
in Ref. 0 as 
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Taking the limit /3 —> 00 , Eq. (fTTH) reduces to 
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In the non-extended phase space, the first law of thermodynamics reads 

dM = TdS + $dQ. 

So the entropy can be derived as 


f“ r + 


S = 


1 (dM ^ ^ Efc(n— l)r" 5 f r + 2kr^_a ^ k 2 a 


T \dr 4 


n—1 n — 3 n — 5 


(16) 


(17) 


(18) 


(19) 


The above result is derived for n > 5 while the integration is divergent for n ^ 5. The charge Q is related to the 
parameter q by 
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Then the expression of the mass can be reorganized as 
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FIG. 1: (a) C<s> vs. r+ for k = 1, n = 6, a = 1, <f> = 1, A = —2 (b) T vs. for k = 1, n = 6, a = 1, <f> = 1, A = —2 


Utilizing Eqs. m and m, the electric potential can be calculated as 
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To study the phase transition in the grand canonical ensemble, it is more convenient to express the mass into the 
function of the electric potential $ as follows 
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The specific heat at constant electric potential can be obtained as 
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One can easily draw the conclusion that the specific heat at constant electric potential may diverge when 


B{r+,*) = 0 , (28) 

implying the existence of phase transition. The above equation can be solved numerically and the results for k = 1,-1 
are presented in Table H1ITT1 respectively, where the impact of the various parameters are studied thoroughly via control 
variate method. It is quite interesting to note that there are two critical points for the case n = 6 , k = 1 while there is 
only one for other cases. And the distance between the two phase transition points becomes larger with the increasing 
of A and 4> while it first becomes larger then becomes smaller with the increasing of a. The case n = 6 , k = 1 is 
shown graphically in FigJTJ Both the behaviors of the specific heat and Hawking temperature are depicted. It is easy 
to find that the two phase transition points where the specific heat diverges are physical for the Hawking temperature 
is positive. The black holes can be divided into three phases. Namely, small stable (C$ > 0) black hole, medium 
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unstable (C$ < 0) black hole and large stable (C$ > 0) black hole. For a more comprehensive picture, we also plot 
the three-dimensional figure for the case n = 6, k = 1 in Fig(S]and for the case n = 6, k = — 1 in FigEI 

The case k = 0 is quite simple. When k = 0, Eq. m can be simplified as 

B(r+,$) = -24Ar®_ + 24(n - 2) 2 3> 2 > 0. (29) 

So there exists no phase transition for k = 0. 


III. THE NATURE OF PHASE TRANSITION IN THE GRAND CANONICAL ENSEMBLE 


In the extended space, it is convenient to utilize the classical Ehrenfest equations to study the nature of phase 
transition at the critical point. However, here, in the non-extended phase space, we would like to introduce the novel 
analog form of Ehrenfest equations proposed by Banerjee et al. [37] as follow 
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FIG. 3: (a) C<s> vs. r+ for k = —1 , n = 6, $ = 1, A = —2 (b) T vs. r+ for k = —1 ,n = 6, a = 1, A = —2 (c) C<s> vs. r+ for 
k = — l,n = 6, a = 1, $ = 1 


TABLE I: The location of phase transition points for k = 1 


n 

a 

$ 

A 

r+i 

r +2 

6 
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1 

-2 

0.406 

1.248 

6 

i 

1 
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0.420 

0.964 

6 

i 

1 
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0.440 

0.794 

6 

0.5 

1 

-2 

0.283 

1.076 

6 

2 
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0.595 

1.363 

6 

1 
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1.592 

6 

1 
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1.666 

7 

1 
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- 

8 

1 
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- 
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1 

1 

-2 

1.812 

- 
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TABLE II: The location of phase transition points for k = — 1 


n 

a 
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1 

1 
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0.386 


a kj 




FIG. 4: (a) a vs. r+ for k = 1, n = 6, a = 1, $ = 1, A = —2 (b) kt vs. r+ for k = 1, n = 6, a = 1, <f> = 1, A = —2 


where a = kt = <j(§$)t are the analog of volume expansion coefficient and isothermal compressibility 

respectively. Their explicit forms can be calculated as follows 
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a, Kj 1 may also diverge at the phase transition point because they share the same factor as C$ in their denominators. 
It can be clearly seen in Fig. 3] 

From the definitions of a and C$, one can obtain 
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So the first equation of Erhenfest equations has been proved to be valid. 

The L.H.S of Eg. (1311) can be obtained as 
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Note that we have utilized the phase transition condition = 0. From the thermodynamic identity [41 
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Noting that in the above derivation, we have also utilized both the definitions of kt and a. We can obtain 
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From Eas. (13711 and fiUl) . one can easily draw the conclusion that the second equation of Ehrenfest equations also 
holds. The Prigogine-Defay(PD) ratio can be calculated as 
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Eg. (TUT) and the validity of Ehrenfest equations show that Lovelock AdS black holes in grand canonical ensemble 
undergo second order phase transition. 


IV. THERMODYNAMIC GEOMETRY OF LOVELOCK ADS BLACK HOLES 


Weinhold’s metric [47| and Ruppeiner’s metric |48j] are defined respectively as 
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And they are conformally connected to each other through the map [63] 
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R 



FIG. 5: Ruppeiner scalar curvature R vs. r+ for k = 1, n = 6, a = 1, $ = 1, A = —2 
Utilizing Eas. (fT71) . (PHI) . (PP5l) . CC>1) and (HT1) . the components of Ruppeiner’s metric can be derived as 
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Utilizing Eas. (l49l) - (l5T1) . we can obtain Ruppeiner scalar curvature as 
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where 
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and E(r+, Q) is too lengthy to be displayed here. The above result has been rewritten in the function of $ so that we 
can compare it with the specific heat. It is not difficult to observe from Eg. (1521) that in the denominator of Ruppeiner 
scalar curvature, the fifth factor is exactly one quarter of the denominator of the specific heat while the last factor 
coincides with the numerator of the Hawking temperature. In other words, the Ruppeiner scalar curvature may 
diverge exactly where the specific heat diverges. It also reveals the extremal black hole condition that the Hawking 
temperature is zero. For an intuitive understanding, one can observe the behavior of Ruppeiner scalar curvature R 
in Figj5] Comparing Fig[5]with Fig 1(a) one can find that the divergence structures of both the Ruppeiner scalar 


curvature and the specific heat are exactly the same. And the Ruppeiner metric does provide a wonderful tool for 
one to probe the phase structures of black holes. 


Among thermodynamic geometry theories, Ruppeiner geometry has been proved to be outstanding for its pro¬ 
found physical meaning. As argued in Ref. [501 ]. Ruppeiner scalar curvature R results from the thermodynamic 
information metric giving thermodynamic fluctuations and may be interpreted physically as the measurement of the 
correlation between fluctuating Planck length pixels of event horizon. In the region with positive R , repulsive interac¬ 
tions (Fermionic behavior) dominate while in the region with negative R , attractive interactions (Bosonic behavior) 
dominate. Moreover, | R \ indicates the average size of fluctuations. 
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V. CONCLUDING REMARKS 


In this paper, we extend our former research of charged topological Lovelock AdS black holes to the non-extended 
phase space. Specifically, we investigate phase transition of (n+l)-dimensional Lovelock AdS black holes in the grand 
canonical ensemble. Firstly, we calculated the specific heat at constant electric potential. To probe the impact 
of the various parameters, we utilize the control variate method and solve the phase transition condition equation 
numerically for the case k = 1, —1. There are two critical points for the case n = 6, k = 1 while there is only one for 
other cases. And the distance between the two phase transition points becomes larger with the increasing of A and 
$ while it first becomes larger then becomes smaller with the increasing of a. We also study the behavior of specific 
heat graphically. As can be seen from the graph, the black holes can be divided into three phases. Namely, small 
stable (C$ > 0) black hole, medium unstable (C$ < 0) black hole and large stable ((7$ > 0) black hole. The graph of 
Hawking temperature is also depicted to check whether the phase transition points locate in the physical region. For 
k = 0, there exists no phase transition point. 

To figure out the nature of phase transition in the grand canonical ensemble, we carry out an analytic check of 
the analog form of Elirenfest equations proposed by Banerjee et al. It is proved that the two Ehrenfest equations 
hold at the phase transition point. Prigogine-Defay ratio is also calculated. Based on these results, one can draw the 
conclusion that Lovelock AdS black holes in grand canonical ensemble undergo a second order phase transition. 

To examine the phase structure in the grand canonical ensemble, we also utilize the thermodynamic geometry 
method. Specifically, we calculate both the Weinhold metric and Ruppeiner metric. It is shown that in the denomina¬ 
tor of Ruppeiner scalar curvature, the fifth factor is exactly one quarter of the denominator of the specific heat while 
the last factor coincides with the numerator of the Hawking temperature. So the Ruppeiner scalar curvature may 
diverge exactly where the specific heat diverges. It also reveals the extremal black hole condition that the Hawking 
temperature is zero. From the graph of Ruppeiner scalar curvature, one can see clearly that the divergence structures 
of both the Ruppeiner scalar curvature and the specific heat are exactly the same. Our research provides one more 
example that Ruppeiner metric serves as a wonderful tool to probe the phase structures of black holes. 

Note that one may vary the spatial dimension, the cosmological constant, and the coefficients of the curvature 
terms in the Lagrangian and we mainly concentrate on a few instances of a very large model in this paper. The 
control variate method has been utilized to crack down the problem of probing the impact of the various parameters. 
We choose such parameter regions that we can compare our results with those in former literatures. One can easily 
extend our results to more cases. The black hole solution here was derived for the special case that the second and 
third order Lovelock coefficients satisfy certain conditions. Phase transition in the non-extended space of more general 
black hole solutions in Lovelock gravtiy would be further investigated in our future work. Also note that the methods 
utilized in this paper can be generalized to an arbitrary non-linear electrodynamics Lagrangian, the specific results in 
this paper however are model-dependent. For a more general analysis, we would like to draw the readers’ attention 
to the excellent work (64J, where the authors presented an elegant procedure for Gauss-Bonnet gravity regardless of 
the explicit form of the nonlinear electrodynamics Lagrangian. It certainly deserves to extend this treatment to the 
third-order Lovelock case in future research. 
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